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Abstract 

^P . The Schwinger-Dyson formalism is used to check the consistency 

of instanton model solutions for the quark propagator with recent 
models of confining gluon propagators. We find that the models are 
not consistent. A major discrepancy is the absence of a vector con- 
densate in the instanton model that is present in the solutions with 
qj ' nonperturbative confining gluons. 
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1 Introduction 

The quark propagator contains valuable information about nonperturbative 
quantum chromodynamics (NPQCD). In the Dyson- Schwinger formalism for 
QCD the integral equation for the quark propagator involves the dressed 
gluon propagator and quark-gluon vertex function, which leads to a set of 
coupled integral equations. Recently this formalism has been used to study 
the quark and gluon propagators. By using theoretical arguments to express 
the quark-gluon vertex in terms of the quark propagator |I], Q or other ap- 
proximations or models for the dressed vertex function it has been possible 
to constrain the nonperturbative gluon propagator by the quark condensates. 
The quark condensates, which arise from chiral symmetry breaking ||, have 
been determined by extensive fits to hadronic properties in the QCD sum 
rule method |4[]. It was recently shown that space-time structure of the con- 
densates || as well as the mixed quark condensate || can be consistently 
determined in the Dyson-Schwinger form with a confining gluon propagator. 
The nonperturbative quark propagator has also been modeled by treating 
quarks propagating in the instanton medium. Instanton models were used to 
attempt to understand the various condensates from the very early days of 
the QCD sum rule method J3J. Based on the known instanton zero modes M, 
an instanton medium picture was shown [|J to be consistent with the known 
quark and gluon condensates. Subsequently, the propagation of light quarks 
in the instanton medium was formulated quantitatively H and refined in its 
treatment of non-zero-mode propagation by Pobylitsa [IIJ. Although it does 
not give rise to a long-range confining force between quarks, the instanton 
vacuum has been shown to provide a good phenomenological description of 
many hadronic properties |T1| , as do QCD sum rules without instantons. 



Instanton effects have been explicitly included in the nucleon correlator in 
QCD sum rule calculations of the nucleon mass \T2, pi and the nucleon 



magnetic moments |l ! 



In the present work we use the forms for the nonperturbative quark propa- 
gator derived in the instanton models in the Dyson-Schwinger formalism and 
check the consistency of explicit gluonic compared to the instanton model 
treatments of NPQCD. We avoid the discussion of including instanton along 
with other nonperturbative QCD effects, such as condensates and Goldstone 
bosons, without double-counting, by simply seeing if the instanton model 
for the quark propagator gives a consistent solution if one uses the dressed 



gluon propagators found in the previous Dyson- Schwinger models within the 
Dyson-Schwinger framework. 

2 Dyson-Schwinger Form for Quark Propa- 
gator 

The quark propagator is defined by 

S q (x) = <0|T[g(x)g(0)]|0>, (1) 

where q(x) is the quark field and T the time-ordering operator. The quark 
self-energy, £, is defined by 

S q (p)- 1 = ip + m q + E(p), (2) 

with m g the current quark mass and the notation ft = J2a1aP a , where the 
7 a are Dirac matrices. 

Starting with the QCD Lagrangian, 

L(x) = qf(x)(iD-m f )q f (x) + L 9lue (x), (3) 

with D Q = d a -i g s A a (x) and A(x) is the QCD color field, one can derive the 



Dyson-Schwinger equation for the quark propagator [|15| . The self-energy of 
the quark is given by 

S(p) = / J^9 2 S D? U (P - Qh^S q (q)l*(q,p), (4) 

with D® b u (q) the gluon propagator, A" the color SU(3) matrix and T^(q,p) the 
quark-gluon vertex. In the present work we use the approximation T^(q,p) = 
"Yu-fi the so-called "rainbow approximation", to be consistent with Ref[Q. It 
is quite straight-forward for us to go beyond the rainbow approximation. 

Eq. (g) expresses the inverse quark propagator as a perturbative and a 
nonperturbative part. Alternatively one can write the propagator as 

S q (x) = S q DT (x) + S^(x). (5) 



For short distances, the operator product expansion for the scalar part of 
S^ p (x) gives 

S? P (x) ~ <: q(x)q(0) :> (6) 

,2 

= <: q(Q)q(Q) :> ~ < 0| : g(0)<x • G(0)g(0) : |0 > + . . . , 

in which the local operators of the expansion are the quark condensate, the 
mixed condensate, and so forth. G c is the field tensor associated with the 
color field. 

We use the Feynman-like gauge and write the gluon propagator as 

DfM = ^^D{q). (7) 

In the present work we use the results of Ref . || for our model of the function 
D(q). 

An important observation is that the inverse quark propagator can be 
written in Euclidean space as 

S q (p)- l = ipA(p 2 ) + B(p 2 ) (8) 

Except for the current quark mass and perturbative corrections, the func- 
tions A(p 2 )-1 and B(p 2 ) are nonperturbative quantities which we refer to 
as the vector and scalar propagator condensates, respectively. The Dyson- 
Schwinger equations (in the Feynman gauge) for A and B are the coupled 

set 

\A( 2a 11 2 8 2 f d*q A(q 2 ) 

B( "» = ¥ s '/(^-« W)+W) - (9) 

The quark and mixed quark condensates are given by 

<: ^ (0):> = -i/^^rW (10) 

9 r J ,B{8){2-A{a)) 



< 0| : q{0)ga ■ G(0)q(0) : |0 > = — f dss[s 



sA 2 (s)+B 2 (s) 
81B(s)[2sA{s)(A(s) - 1) + B 2 {s)}. 
16{sA 2 {s) + B 2 {s)) 



These equations are all obtained in the rainbow approximation. The expres- 
sion for the mixed quark condensate was derived in a 1/N C expansion of the 
gluon two-point function||, making use of the self-consistency relations (|5|). 
In Ref. the effective quark mass had reached the current quark mass 
by 1 GeV, which was taken as the renormalizaton point and cutoff for the 
condensate integrals. We do not use such a cutoff, consistent with Ref. [|I] 



We use the notation &=-(2ir) <: g(0)g(0) :> for the quark condensate and 
m 2 a = (27r) 2 < 0| : q(0)gcr ■ G(0)g(0) : |0 > for the mixed quark condensate. 
The values of these condensates have been determined by fits to experiment 
in QCD sum rule calculations to be a ~ 0.55 GeV 3 and m 2 ~ 0.8 GeV 2 , or 
m 2 a ~ 0.44 GeV 5 . 

The coupled integral equations Eq.(^) have been solved for with various 
forms and the parameters fit to the quark condensate. It was shown in Ref || 
that this gave a satisfactory result for the mixed quark condensate and in 
Ref|| that the nonlocal scale of the nonperturbative quark propagator as 
determined from fits to deep inelastic scattering [pj]] was also fit satisfactorily. 
The range of the confining infra-red gluon propagator, which is taken as a 
Gaussian in this work, is very short. In other calculations |L6| the confining 
part of the gluon propagator had a longer range, but the fit to the quark 
condensate was not very good. We also use these models for D(s) in our 
present study. 

3 Instanton Model for Quark Propagator 



In the instanton calculation of Ref |T0[ the inverse propagator is of the form 
(for the current quark mass m q = 0) 

S^p)' 1 = tp + Bjip 2 ) 
Bj(p) = Bl + OiN/VNc). (11) 

Using the leading term of a series expansion, the author found the following 
instanton solution for the inverse quark propagator, Eq.(|8|): 

A I {p) - 1.0 = 0.0 

Bj(p) = Kp 2 f{\p) 
o 

f(p) = --(3I (p) + I 2 (p))xK 1 (p), (12) 



where K = 0.29 GeV -1 , and p is in units of GeV. We refer to Eq. (Jf^) as 
the Pobylitsa solution. 

An important feature is that the vector propagator condensate vanishes: 
A(s) - 1 = 0. This follows from the symmetries of the model used by Pobyl- 
itsa. 

This solution gives an effective quark mass which falls as 1/p 6 for p larger 
than a few hundred MeV. Other treatments of the effective quark mass show 
different behavior. For example a light-cone model for the pion form factor 
[ PI gives B(s) as almost constant and equal to the constituent quark mass 
for p less that about 1 GeV and dropping rapidly to zero at about 1 GeV. 
In this model A ~ 1.0. Let us consider such a model, which we call the Ic 
model: 

A lc (p)- 1.0 = 0.0 

B lc = ^ for p 2 < So (13) 

= 0.0 for p 2 > s , 

where s is a constant to be determined. This Ic model can be considered 
as a simple alternative to the Pobylitsa model in the sense that the vector 
propagator condensate vanishes. 

The mixed quark condensate for the instanton model of Ref || was 
estimated |19| to be m 2 = 1.4 GeV 2 , about a factor of two larger than the 
empirical QCD sum rule value. 

4 Results and discussion 



We use for the form of the gluon propagator ||16[ D(s) 

„ . . 9 V 2 s An 2 d . . 

* D M = 3,r V~ I + Stated" (14) 

with the parameter c„ = (1.0,0.0) to (include, neglect) the perturbative ultra- 
violet behavior. The strength parameter X and range parameter A are deter- 
mined by solving the coupled Dyson- Schwinger equations, Eq.@, and fitting 
f n , the pion decay constant, and the quark condensate through Eq. (fl0|) . For 
the Feynman gauge with c u = 0.0 these parameters are |X = 1.4 GeV and 
A = 2.0 x 10~ 3 GeV 2 . 

6 



For the Ic model, using a = -(27r) 2 <: q(0)q(0) :> = .55 GeV 3 , we find s 
= 0.8 GeV 2 

m 2 a = 0.831 GeV 5 
B(0) = 4.64 GeV 
A(0) - 1.0 = 7.28 (15) 

This solution clearly gives inconsistent values for A and B and the mixed 
gluon condensate is about a factor of 2 larger than that generally found in 
QCD sum rules. Note that this value for the mixed condensate is about the 



same as that obtained[|T9j with an instanton model |J. 

For the instanton model, Eq. flT2|) , the value of K found by Pobylitsa gives 
the value of B/(0) ~ 420 MeV, rather than the constituent quark mass. By 



modifying the factor of K in Eq. (|l2[) by a factor of Mtv/(3 x .420) so that 
B/(0) = Mjv/3, we obtain a value of quark condensate in agreement with the 
phenomenological one, and find for the instanton model 



a - 


= 0.60GW 3 


m 2 a - 


= 2.04 GeV 5 


B(Q) - 


= 4.74 GeV 


A(0)-1.0 = 


= 7.71. 



(16) 

Therefore we find that the Ic and Pobylitsa solutions give similar results 
for the Dyson-Schwinger formalism for the quark propagator. Although the 
value for the quark condensate agrees with the empirical one, the mixed 
condensate is almost a factor of 5 larger than the empirical value of .44 
GeV 5 , and about 2.6 times larger than the estimate of Ref. [jT9| . 

The range of the infra-red part of the gluon propagator is very small 
in these calculations in which the gluon propagator is constrained by the 
phenomenological value of the quark condensate, a=.55 GeV 3 . If we relax 



this constraint we can use the models of Ref. |16j in which longer-range 
infra-red parts of D(s) were used. The results for the Pobylitsa model for 
X=1.53 GeV and A = 0.02 GeV 2 with c u = 1.0 and A = 0.2 GeV are 

B(0) = 5.92 GeV 
A(0)-1.0 = 9.89; (17) 



and with X = 1.65 GeV and A = 0.2 GeV 2 

B(0) = 1.83 GeV 
A(Q) - 1.0 = 3.14. (18) 

Recall that for all of these quark propagator solutions the value of B(0) is 
the value of the constituent quark mass, taken as .313 GeV. From Eq . (|18|) 
we note that for the longest range gluonic confining propagator, with a range 
of about 0.8 GeV, the instanton solution is most nearly consistent, but the 



value of the quark condensate is quite small in that model [16 . 

Our results show that the instanton model is not consistent with the 
Dyson-Schwinger model based on the rainbow approximation using a con- 
fining gluon propagator chosen to fit the condensates and other hadronic 
properties. One major difference between the NPQCD calculation solving 
the model with a confining gluon propagator and the simple substitution of 



the instanton model of the quark propagator Ref|10|| in the Dyson-Schwinger 
integrals is the vanishing of the vector propagator condensate in the instan- 
ton model. In future work we plan to investigate the inclusion of instanton 
effects along with modified confining gluon propagators. 

This work was supported in part by the National Science Foundation 
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